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ABSTRACT 

LS 5039 is a high-mass binary with a period of 4 days, containing a compact object 
and an O star, one of the few high-mass binaries detected in 7-rays. Our Chandra ACIS 
observation of LS 5039 provided a high-significance (~ 10cr) detection of extended emis- 
sion clearly visible for up to V from the point source. The spectrum of this emission can 
be described by an absorbed power-law model with photon index V = 1.9 ± 0.3, some- 
what softer than the point source spectrum V = 1.44 ± 0.07, with the same absorption, 
Njj = (6.4 ± 0.6) x 10 21 cm~ 2 . The observed 0.5-8 keV flux of the extended emission is 
~ 8.8 x 10~ 14 ergs _1 cm -2 , or 5% of the point source flux; the latter is a factor of ~ 2 
lower than the lowest flux detected so far. Fainter extended emission with comparable 
flux and a softer (r ~ 3) spectrum is detected at even greater radii (up to 2'). Two 
possible interpretations of the extended emission are a dust scattering halo and a syn- 
chrotron nebula powered by energetic particles escaping the binary. We discuss both of 
these scenarios and favor the nebula interpretation, although some dust contribution is 
possible. We have also found transient sources located within a narrow stripe south of 
LS 5039. We discuss the likelihood of these sources to be related to LS 5039. 

Subject headings: binaries: individual (LS 5039) — stars: winds, outflows — X-rays: 
binaries 



1. Introduction 

In a high-mass X-ray binary (HMXB), a compact object, which can be a black hole (BH) or a 
neutron star (NS), accretes some portion of the wind from a high-mass (M > 5Mq) non-degenerate 
companion star. Some HMXBs are called 'micro-quasars' (/iQSOs) when resolved, collimated radio 
features are observed, which are interpreted as relativistic jets. 
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There has been renewed intere st in /uQSOs since t he detection of parsec-scale X-ray jets ([Corbel 



20071 ) and variable 7-ray emission (jAbdo et al.ll2009bl ). The nature of the high-energy emission for 
two micro-quasars, LS 5039 and L S I +61 303, has been m uch debated, e.g., the emission could 
truly be generated in a /xQSO jet (jKhangulvan et al.l 120081 ) . or possibl y it could be generated by 



Dubuil2Q0J). 



the interaction of a young pulsar wind with the wind of the companion (|Leahvll2004l ; iDhawan et al 
200fil: 



LS 5039 was first identified as an HMXB by iMotch et al.l (|l997l ) based on the positional co- 
incidence of an X-ray source with an 07V star (V = 11.5), proje cted near the Gala ctic bulge 
(7 = 16?8816, b = — 1?2892). Using optical and infra-red spectroscopy. IClark et al.1 (120011) found the 
companion star to be type 06.5V(f), variable in the IR but not in the optical. iMcSwain et al.l (|200ll ) 
measured line radial velocities and found binary period P or b = 4. 1 d and substantial eccentricity 
e = OA. Using deta iled time-resolved spectroscopy and atmospheric modeling of the companion, 
Casares et al.l (|2005l ) refined the binary parameters (-P rb = 3.906 d, e = 0.35) and concluded that 
the system must harbor a high-mass compact object, probabl y a BH , for a n inferred distance of 
2.5 ± 0.1 kpc. The conclusion was, however, later disputed by iDubusI (|2006l ) who argued that the 
NS option still cannot be dismissed. Using high-pr ecision, long time base-line photometry from 
the MOST satellite and ground-based spectroscopy, ISartv et al.l (|201ll ) further refined the orbital 
characteristics, preferring a slightly lower eccentricity e = 0.24 ±0.08 and component masses 26M 
for the O-star and 1.7M Q for the compact object. Thus, the nature of the compact object in LS 
5039 remains unknown. 



An unresolved, non-thermal, variable radio source was fist detected by iMarti et al.l (|l998l ). 
Using the high ang ular resolution of VLB A, the source was resol ved, and 'jets' o bserved by 
Paredes et al.l (|2000l ); the source was therefore classified as a /uQSO. iRibo et al.l ([2002]) measured 
the system proper motion (~12masyr _1 , using radio and optical astrometry) and, based on its 
high tangential velocity, called the system 'runaway', i.e., with enough velocity to leave the Galac- 
tic disk. The proper motion vecto r can be traced ba ck to the supernova remnant (SNR) G16.8— 1.1. 
Subsequent radio observations by lRibo et al.l (|2008l ) showed that the resolved emission of LS 5039 
varies in intensity, and that the position angle of the elongated features changes. 

LS 5039 has been extensively observed in X-rays, mainly to study its variability and spectrum. 



The X-ray emission is modulated at the orbital period. The light-curve 



data taken in 1999-2007 was found to be remarkably stable (Kishishit a et al 



aperiodic variability has been reported earlier with RXTE (|Reig et al 



2003 



produ ced from X-ray 



2009 ), although some 



Ribo et al.lll999l ). 



Both LS 5039 and LS I +61 303 have been detected at higher energies. L S 5039 has been de- 
tected with HESS (lAharonian et al.lhoOfih . INTEGRAL (|Hoffmann et alJbood ) and Fermi (|Abdo et al 
2009al ). with a spectrum similar to 7-ray pulsars. The hi gh-energy light curves show modula tions 
at the orbital period, peaking at superior conjunction (see iBosch- Ramon Khangulvanll2009l for a 



review . 



In this paper we report the results of a 38 ks imaging observation of LS 5039 with Chandra. 
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In Section [2j we describe our data and reduction procedures, as well as several archival data 
sets. In Section [3] we analyze the spectra of both the point source and the extended emission, 
calculate radial profiles, and examine the extended morphology. We also describe the properties 
of a number of transient sources that appear in the field. In Section H] we discuss our findings and 
their implications. We conclude with a brief summary in Section [5j We also include Appendix A, 
in which we calculate the X-ray contribution due to dust scattering in the direction of LS 5039. 



2. Observations and data analysis 

Table [JJ lists the X-ray observations considered in this work, including five archival observations 
(from XMM-Newton and Chandra) and our latest Chandra observation. 



2.1. New Chandra observation 



LS 5039 was observed with the Advanced CCD Imaging Spectrometer (ACIS) on board Chan- 
dra on 2009 October 30 (ObsID 10696). Ou r observation spanned orbital phases cf> = 0.064-0.184 
(using the ephemeris of ICasares et alJ 120051 . where 4> = corresponds to periastron). The target 
was observed for 38.55 ks in timed exposure mode using "very faint" telemetry format. To minimize 
pile- up effects, we used a custom 750 pixel wide sub- array and turned off all the ACIS chips ex- 
cept 13. The frame time in this configuration is 2.24104 s (2.2 s exposure time and 0.04104 s image 
transfer time). The Y-offset was — 2'.8, placing the source 3'.03 off axis. This resulted in a broader 
PSF, further helping to mitigate pile-up. There were no significant background flares during the 
observation. The useful scientific exposure time (live-time) was 37.84 ks. 

The data were reduced and analyzed with the Chandra Interactive Analysis Observations 
(CIAO) package (ver. 4.2), using CALDB 4.2.0. To produce images of the pulsar and its vicinity at 
subpixel resolution, we removed the pip eline pixel random ization and applied the subpixel resolution 
tool to split-pixel events in the image ([Mori et al.ll200ll ). 



The image produced from this data is shown as the bottom-right panel in Figure [TJ 



2.2. Archival observations 



We downloaded archival observations of LS 5039, three by Chandra and two by XMM-Newton, 
see Table [JJ In the case of the Chandra data, we used the standard Level 2 calibrated science event 
files. We used the energy range 0.59-6.8 keV, which is the optimal range for viewing the extended 
emission (see below) , and show the resultant images in Figure [TJ Although taken with the same 
instrument, the data sets span a range of exposure times and off-axis distances. 
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Fig. 1. — Images of the field of LS5039. All images are 7'x7', in the 0.59-6.8 keV energy band for 
Chandra, 1-7 keV for XMM-Newton. In the bottom-right panel we mark the source and background 
extraction regions used in our spectral analysis (green circles; see §3.2). White numbered circles 
show the faint field point sources discussed in §3.31 
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In the case of the XMM-Newton data, we made use of only the EPIC-PN data, from the 
standard pipeline in the XMM-SAS calibration package. These data were filtered in the energy 
range 1-7 keV, and the resultant images are also shown in Figure [TJ 

We performed no further processing on any of these archival data, and use them only for 
comparison purposes while looking for transient sources ( §3,3p . 

In addition to the six observations discussed in this paper, there were two observations by 
Chandra and two by XMM-Newton. The Chandra data were taken in continuous-clocking mode, 
and so are not useful for analysis either of extended emission or faint field sources because of the 
much higher background. For the XMM observations, the spacecraft orientation was such that the 
locations of faint field sources shown in Figure [1] were not imaged, and again the PSF is too broad 
and the background too high to see extended emission. We consider neither observation further, 
except to note the point source flux measured by Chandra, below. 



From the images in Figure HJ it appears that LS 5039 is a bright unresolved point source with 
extended emission around it; the extended emission is seen most clearly in the latest Chandra 
observation (bottom right panel). We split our results into three sections: the spectrum of the 
bright point source, which is essentially unaffected by faint extended emission, the properties of the 
extended emission, which includes some contribution by the bright source's PSF, and a description 
of the transient sources in the field. 



From the Chandra ACIS data of 2009 October 30, we extracted the spectrum of the point-like 
source corresponding to LS 5039, centered at coordinate f] R.A. = 276?56278, deck = -14?84871 
(J2000), using the CIAO task psextract. We picked background regions far from the source, and 
as free of point sources as possible, yet large enough (total area 51,100 arcsec 2 ) that the background 
spectrum be well-determined. The extraction regions are shown on Figure [H 

We fitted an absorbed power-law (PL) model to the spectrum, extracted from the r = 5" 
aperture, using the Sherpa package. We grouped the counts in 25-count bins and fitted in the 
energy range 0.5-8 keV (4602 counts in total, of which 3.8 counts belong to the background; source 
count rate 0.122 ± 0.002 counts s _1 ). The observed energy flux is F 5 _ 8]ie y = (1.92 ± 0.03) x 
10~ 12 erg cm~ 2 s^ 1 . The fitting parameters are given in Table 2, and the fluxed count-rate spectrunj^ 




3. 



Results 



3.1. 



Point source 




x We use the coordinate system provided by the data processing pipeline. 

2 The fluxed count-rate spectrum, which we will call the 'fluxed spectrum' for brevity, is calculated by dividing the 
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is shown in Figure [2j 



The inferred extinction-corrected flux for the 1-10 keV range is 2.6x10 12 erg cm 2 s . 



Al- 



though the phase interval of our observation corresponds to the minimum of the previously observed 
orbital variability, t his is a factor of two lower than any flux that h as been noted in the past (e.g., the 
long-term study by iKishishita et al.ll2009l : iTakahashi et al.ll2009l ). For two Chandra observations a 
few months before, 2009 July 31 (observation ID 10053) and August 6 (observation ID 10932; PI 
NandaRea), with orbital phases <f) =0.81-0.95 and 0.35-0.43, we found fluxes 8. 7x 10 -12 ergcm _2 s _1 
and 7.7xl0 -12 erg cm -2 s _1 , respectively, close to the values previously measured at these phases. 

With a count rate of 0.27 counts per frame, the pile-up fraction would be 10% if the target 
were on-axis (estimated using Chandra PIMMS). As the source was ~ 3' off-axis in our observation, 
the actual pile-up fraction is about 4%, as we estimated with the aid of MARX simulations (see 
below). 

Our best-fit photon index, T = 1.44 ±0.07, is only marginally consis tent with those previously 
measured at the same orbital phase, V =1.55-1.6 (IKishishita et al.ll2009i ). 



3.2. Extended emission 

We used Chandra Ray Tracer (ChaRTj! and MARX0 software to simulate the PSF for the 
given spectrum and location on the detector. In order to reduce the Poissonian noise, we made 
a simulated event file for 200 ks exposure, filtered the same way as for real data, and then scaled 
by the exposure times. We used ChaRT with the input spectrum derived from our spectral fit 
to the central point source and MARX Version 4.4 with parameters ACIS_Exposure_Time=2 . 2 (to 
account for the non-standard sub-array mode) and DitherBlur=0 . 27 (a measure of the aspect 
reconstruction accuracy and pixelization by the detector, typically near / /3 for ACIS). All other 
parameters were left at their default values. 

The simulated PSF is shown in Figure HI We calculated the signal-to-noise ratio (S/N) of the 
extended emission near the point source, and found the optimum energy range 0.59-6.8 keV. In the 
radii range 20"-60", we find 745 counts, of which expect 378 from the background and 81 from 
the PSF wings, giving 286 ± 29 source counts, a detection significance of 9.9cr. In the radii range 
60"-120", there are 1766 counts, of which 1278 should be background and 32 from the PSF wings, 
giving 456 ± 51 source counts, a detection significance of about 8.9<r. 



count rate in a given energy bin by an average effective area corresponding to this bin. 
3 See http : //cxc .harvard. edu/chart /threads/ index.html 
4 See |http : //space .mit . edu/CXC/MARX/ 
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Table 1. Log of observations 



Date 


ObsID 


Telescope/ 


Off-axis 


Exposure 






instrument 


angle (') 


time (ks) 


2004-07-09 


4600 


Chandra/ACIS 


12 


11.0 


2004-07-11 


5341 


Chandra/ACIS 


12 


18.0 


2005-04-13 


6259 


Chandra/ACIS 





5.0 


2005-09-22 


0202950201 


XMM/EPIC 





15.8 


2005-09-24 


0202950301 


XMM/EPIC 





10.4 


2009-10-30 


10696 


Chandra/ACIS 


3 


37.8 



Note. — Exposure times are live (i.e., corrected for GTI filtering and 
dead-time) . 



Table 2. Power-law spectral fits 



Region 


Radius range (") 


N H (cm" 2 ) 


r 


A a 


X 2 /dof 


punabs b 
^0.5-8keV 


Point source 


0-5 


6.4(6) x fO 21 


1.44(7) 


35.1(3) 


139.3/176 


25.4(2) 


Extended 


20-60 


6.4xl0 21 c 


1.9(7) 


2.9(7) 


47/40 


1.4(3) 


Extended 


60-f20 


6.4xl0 21 c 


3.1(5) 


8(2) 


70/103 


2.4(6) 



Note. — Fits were performed in the 0.5-8 keV energy range. Numbers in parentheses are la 
uncertainties on the last digit. Confidence contours are shown in Figure [3l 

a Normalization in 10~ 5 photons cm~ 2 s _1 keV _1 at IkeV. 

b Unabsorbed energy flux in 10~ 13 ergcm~ 2 s -1 . 

c Held fixed in the fit. 
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Fig. 2. — Fluxed spectrum of LS 5039 (black) and the extended emission extracted from 20"-60" 
(red) and 60"— 120" (green) annuli around LS 5039. Best-fit absorbed power-law models are over- 
plotted as solid lines, and residuals, x = (data — model)/uncertainty are shown in the lower three 
panels. The best-fit model and fitting parameters were estimated using standard forward-fitting 
(i.e., by fitting the count-rate spectrum). 
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Fig. 3. — Confidence contours at one-, two- and three-sigma, for the PL fit to the spectra. The 

fits to the annular spectra have Nh fixed to the value obtained for the point source. Lines of 

constant unabsorbed flux in the 0.5-8 keV band are shown with dotted lines, and labeled in units 
er g g -i cm - 
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Fig. 4. — Top: Radial intensity plot around the point source, for the total (gross), background 
subtracted (net), simulated PSF count intensities and dust scattering halo model in the energy 
range 0.59-6.8 keV. Bottom: The difference between the net and PSF curves, i.e., the excess emission 
around the point source. 



-11- 




-i 0.045 

- 0.040 

- 0.035 

- 0.030 
a 0.025 
1 0.020 

0.015 
0.010 
0.005 




Fig. 5. — 256" x 256" image of the vicinity of LS 5039 after subtraction of the simulated PSF. The 
color-bar scale is counts per / .'49 x / .'49 pixel. The image has been smoothed by a Gaussian kernel 
with r = 3". The inset shows the central 27"x27" section, where binning was / /25, and smoothing 
width r = O'.'S. It better demonstrates the residual structure around the point source, likely caused 
by imperfections in the PSF model. 
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Fig. 6. — Histogram of total counts in the r = 20"-60"annulus and optimal energy range (0.59- 
6.8 keV), as a function of azimuthal angle about the point source, with zero angle west and increasing 
towards the north. 
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To analyze the spectrum of the extended emission, we extracted counts from the two annular 
regions, r = 20"-60"and r = 1'— 2', using the CIAO task specextract. The background regions 
were the same as for the point source above (see Figure [1]) . The spectra were once more grouped 
to at least 25 counts per bin, and fitted with an absorbed PL model in the energy range 0.5-8 keV. 
The column density was kept fixed at the best-fit value obtained from the fit to the point source 
spectrum (Njj = 6.4 x 10 21 cm -2 ). For the inner annulus, there is considerable contamination 
from the PSF wings, which may have a different spectrum from the point source because of the 
energy dependence of the the PSF. To account for this, we fit the counts extracted from the 
simulated PSF events in the same annulus, and then include this as an additional frozen component 
(r = 0.59, N = 3.83 x 10~ 7 pns~ 1 cm _2 keV _1 is 1 keV) when fitting the extended emission. The 
results of the fits are shown in Table [21 while the corresponding fluxed spectra and confidence 
contours are shown in Figures [2] and [3l respectively. The observed energy flux was ^0.5-8 keV = 
(8.8 ± 1.9) x 10~ 14 ergcm _2 s _1 (inner) and (7.5 ± 1.9) x 10~ 14 ergcm~ 2 s _1 (outer). Using the 
'optimal' energy range as defined above gives very similar fit parameters. We find that the inner 
annulus spectrum is only slightly softer than the point source, while the outer annulus is significantly 
softer. 

To better see the morphology of the extended emission, we subtracted the model PSF from 
the original image. The relative positions of the simulated and observed sources were adjusted to 
minimize the residuals. The result is shown in Figure [5j Some extended structure is apparent, 
such as wisps to the north and east. With the relatively low S/N it is hard to say more, but, 
interestingly, similar structure is visible in the short on-axis Chandra observation of 2005 April 13 
(Figure [T]) . To quantify any anisotropy, we considered the distribution of the extended emission 
over azimuthal angles about the point source. Figure [6] shows a net histogram of background- and 
PSF-subtracted counts for the inner annulus. The distribution appears to be non-uniform, with 
X 2 = 15.9 for 5 degrees of freedom for a constant value (null hypothesis probability of 0.7%). 

The difference between the simulated and point source images shows residuals up to r = 3", 
although the two agree well in the radial plot (Figure [4]) . This suggests some imperfections in the 
PSF model on small scales. 

3.3. Transient field sources 

A number of small brightness enhancements are apparent in the images in Figure [TJ We have 
labeled the significant ones with white markers. They are all consistent with being point-like, 
although with the small numbers of counts, this is not very restrictive. 

The numbered sources in Figure Q] are projected onto the same sky coordinates in Figure [71 
They appear to lie within a narrow stripe almost due south from LS 5039 in two clumps or clusters. 
All the sources have rather hard spectra (mean energies 2-4 keV). Notably, the source that was 
bright on 2004 July 9 (Source 1) had completely disappeared two days later (even though the latter 
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Fig. 7. — Detected sources in th e field of LS 5039. We show LS5039's proper motion vector (1000 
times yearly value measured by IRibo et al.l |2002| . black arrow) an d the position angles of radio 
elongation at two different epochs (dashed lines; IRibo et al.l 120081 ). The object numbers (bold) 
refer to the white labels in Figure [TJ 
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observation had a longer exposure time, and would have detected it); similarly Source 5 appeared 
on a similar time-scale. We cannot exclude that sources 1 and 5 are the same object, see below. 
The fluxes of these sources are given in Table 

We show zoomed-in views of the locations of the field sources detected during observations 
of LS 5039. In Figure [8] the positions are plot ted, after bore-sight correction, along with sources 
from the 2MASS catalog (jSkrutskie et al.ll2006l ). The nearest 2MASS sources, which have positions 
consistent with (some of) the field sources, are listed in Table H] their infra-red magnitudes. Where 
possible, we also include optical magnitudes. 

Sources 2 and 3 are only marginal detections. If they are not real sources, or if they are 
unrelated to the rest of the sources detected, then it is likely that sources 1 and 5 are a single 
flaring source possibly associated with 2MASS source C. Likewise, sources 4 and 6 could be a single 
persistent source, perhaps associated with 2MASS source A. 
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Fig. 8. — Close-up of the locations of field sources shown in Figure 7, northern group (left) and 
southern group (right). Also shown are sources from the 2MASS catalog (squares). The nearest 
stars are labeled by letters, and their photometric data and positions are given in Table [U The 
coordinates shown here are after the application of a bore-sight correction. Scale bars are 5" in 
length. 



4. Discussion and conclusions 

Our analysis has revealed a highly significant extended emission discernible up to ~ 2' from LS 
5039 (limited by the FOV boundaries). This emission could be either due to the dust scattering of 
X-rays emitted by the unresolved source or it could be emission from relativistic particles escaping 
from the vicinity of LS 5039. Below we discuss both options. In addition, we observed at least 
two faint unresolved sources. These could be field stars, background highly-absorbed CVs or 
quiescent XRBs (or even AGNs seen through the Galactic plane), but they might in principle be 
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Table 3. Fluxes of transient field source fluxes 



Date 


SonrppK 


Flux a 






North 


South 


2004-07-09 


1 


3.6 ±0.8 


<4.4 


2004-07-11 




<0.8 


< 1.4 


2005-04-13 


2,3 


2.4 ±1.1 


2.3 ± 1.7 


2005-09-22 


4 


<0.7 


1.3 ±0.3 


2005-09-24 


5 


3.1 ±0.5 


<2.0 


2009-10-30 


6 


<0.47 


1.6 ±0.4 


a Observed, 




background-subtracted 


energy flux 


(1.0-6. 


8 keV) in 


units of 


10~ 14 erg cm - " 




Uncertainties are at 


1-cr limits are 


at 3-cr. 


Columns 


'North' and 



'South' represent sources in the two apparent 
spatial groupings, see Figures [7] and [SJ 



Table 4. Photometry and positions of 2MASS sources near faint X-ray sources 



Label 


2MASS ID 


J 


H 


K s 


Other 


R.A. (°) 


decl (°) 












photometry 3. 






A 


18261441-1453235 


14.40 


12.62 


11.91 




276.560075 


-14.889871 


B 


18261416-1453275 


13.68 


12.64 


12.22 




276.559035 


-14.890979 


C 


18261389-1452239 


14.40 


>13.50 


>13.28 


B = 16.64, V = 16.20 


276.557879 


-14.873326 


D 


18261419-1452324 


12.50 


11.03 


10.46 


R = 19.4 


276.559134 


-14.875684 



Note. — Labels refer to the those shown in Figured! 2MASS photometric uncertainties are 0.04-0.06 mag. 
a Magnitudes noted in the NOMAD catalog. 
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manifestations of a collimated outflow from LS 5039. 



The flux of the point source LS 5039 was a factor of ~2 lower than expected for its orbita l 
phase, despite previously appearing to have a stable orbital li g ht-curve (IKishishita et al.ll2009l ). 
Previous studies did suggest larger variability (IReig et al.1 120031 : iRibo et al.1 Il999l ) , but the use of 
different instrumentation (particularly RXTE) made the variability hard to verify. A decrease in 
flux could perhaps be explained by a change in accretion state (if the compact source is a BH) but 
even for pulsar emission, it can be explained by changes in the wind pro perties of the com panion 
O-star which also manifest themselves via the variability of H Q emission (|Reig et al.l l2003). 



4.1. Dust scattering halo 



Scattering of X-rays emitted by the unresolved source off the dust lying along the line o f sight 
produces extended emission known as a dust scattering halo (e.g., iPredehl fc Schmittlll995l ). The 
halo brightness scales linearly with the intervening dust column, which is usually assumed to be 
proportional to the hydrogen column density, Nh, measured from X-ray spectra, i.e., scattering 
optical depth at IkeV is r S) ik e v — S(Nh /10 22 cm -2 ), where the proportionality constant S can 
only be determined empirically and depends on propert i es of t he intervening dust which may differ 
for different objects. For instance, IPredehl fc Schmittl (|1995l ) found a mean S ~ 0.5 from their 
sample of X-ray halos, which corresponds r s u^y — 0.24 for LS 5039 where Nh = 6.4 x 10 21 cm -2 
(Table 2). An even better correlation for r s was obtained with the optical extinction t owards the 



Riboetal 



objec t, r Si ikeV = 0.056Ay (see Appendix A), which gives r S) ik e v = 0-22 (for Ay = 3.9 
2002 ) . in good agreement with the previous estimate. The relatively low value of r Sj ik e v allows us 
to neglect the contributions of multiple scatterings for E > 1 keV. 

We have calculated the ha lo profile in the sing le-scattering approximation by folding the spec- 



tral intensity of a halo (see e.g-. lMathis fc Ledll99ll . also Appendix A) with the detector response in 
the 0.5-8 keV energy range. We find that, for instance, for the parameters = 360" and S = 1, and 
the dust distribution function f(x) defined in Appendix A, the dust halo model overall describes 
the observed radial profile (see Fig. [4]). Although some deviations are noticeable, they might be 
attributed to the imperfec t choice of the above parameters (we did not perform rigorous fitting) 
or to the inaccuracy of the iDraind (120031 ) model or the Rayleigh-Gans (RG) approximation to the 
scattering cross section at low energies (< 2 keV). 

The dust distribution function f(x) used to calculate the halo model (see Equation (A8) and 
Figure U]) assumes a lack of absorption/scattering in the immediate vicinity of the source. Such 
dust distribution is consistent with the lack of evidence of significant intrinsic absorption in the 
LS 5039 binary. The upper limit on the intrins ic column density is as low as 2.6 x 10 21 cm -2 
(jBosch-Ramon et al.l 120071 : iTakahashi et al.l 12009 ), while the tota l abso rbing column density mea- 
sured from X-ray spectra changes by < 10% (ITakahashi et al.l l2009l ) with binary phase. Low 
intrinsic absorption would be expected if most of X-rays are emitted far from the donor star, where 
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its wind is sufficiently rarefied (e.g. Szostek k, Dubuslhoill ). 



We should also point out that, in addition to the dust distribution profile, the dust scattering 
model we use has two free parameters (S and 0) which can attain values within rather broad 
ranges, depending on the actual properties of the dust grains. Therefore, the mere fact that the 
model qualitatively fits the observed brightness distribution does not guarantee that the extended 
emission is a dust halo. In fact, the observed azimuthal asymmetry and the hard spectrum of the 
inner extended emission would be difficult to explain by a dust halo. Significant spectral softening is 
expected from halo models. Specifically, for the dust model used, the spatially integrated (between 
r = 20" and 60") spectrum is expected to have r\ a i ^ r po i nt + 2 ~ 3.5, because the scattering cross- 
section a s (E) oc E~ 2 (see Appendix A). Instead, we find a significantly harder spectrum, T ~ 1.9. 
On the other hand, the outer extended emission is indeed much softer and more symmetric, hence 
could have some contribution from a halo. 



4.2. Extended nebula of LS 5039 



The nature of the compact object (BH or NS) in LS 5039 remains subject of a debate. LS 5039 
(and LS I +61 303) are quite different in their temporal and spectral properties from other //QSOs 
and HMXBs, which show transitions between different states and large variations in luminosity. 
These two systems are also, to date, the only //QSOs whose VHE emission has been firmly detected 5 * 



The X-ray light curve of LS5039, which shows remarkable long-term stability (jKishishita et al 



2009 



Takahashi et al. [Hoffmann et al.ll2009 ). and peculiar, variable radio morphology (|Rib6 et al 



20081 ) argue against the accretion scenario. 



If, despite the lack of the usual manifestations of ac cretion, the compact object in LS 5039 is 
a BH accreting in an unusual regime (jCasares et al.ll2005l ). it may still be possible that it produces 
relativistic particles, e.g., via the Blandford-Znajek process (jBlandford fe Znaieklll977l ) or in an 
MHD jet. In this case, an extended nebula could still be formed. For the case of a black hole, 
the majority of hig h-energy emission is expected to be produced in a jet-type axial outflow (e.g., 
Paredes et al.ll2006l ). The p oss ible observational manife stations of such outflows and a discussion of 



their e nergetics is discussed by lRussell &: Fenderl (120101 ) and their extended emission by lHeinz et al 



(|2008l ). The relativistic electrons responsible for the 7-ray emission in the inner jet will cool with 
distance. Flow speeds are typically close to c (lYadavll2006l ). and so X-ray emission ca n be expected 
out to appreciable distances from the black hole, e.g., the plasma clumps seen by ICorbel et al. 
(|2002l ) up to ~ 30" from fiQSO XTE J1550-564. Furthermore, the spin axis of the black hole may 



precess, and so the orientation of the outflow change with time, as with SS 334 (jBlundell et al.ll2007l ) 
and thus fill an extended volume with energetic emission, rather than the more typical collimated 



5 The only other HMXB firmly detected at VHE is B1259— 63, where a young pulsar orbits a Be star on a 3.4 yr, 
highly eccentric orbit. 
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outflows. Such large-angle precession may be hard to explain, however, possibly requiring a third 
massive component to the system. 

A possible alternative interpretation of the high-energy emission from LS 5039 is a pulsar wind 
nebula (PWN) powered by a young, energetic pulsar. In this scenario, modeling of the intra-binary 
shock suggests that the O-star wind dominates dynamically over the pulsar wind, which makes 
the shocked wind flow away from the head of the intra-binary bow-shock and along the bow-shoc k 



surface which is curved away from the O-star (see e.g., Figures 1 and 2 in lSzostek &: Dubusll2011 



If the opening angle of the bow-shock is small, it would resemble a cometary tail rotating around 
the O-star together with the pulsar and pointing in t he direction of the O -star wind (cf. variable 
ra dio morphology of L S I +61 303 resolved with VLBI; iDhawan et abj|2006] ; see also the discussion 



m 



Pavlov et al 



2011a 



Emission from the innermost regions of the bow-shock may be responsible 
for the unresolved, bright X-ray source whose surface brightness varies with the orbital phase 
owing to Doppler boosting (jDubus et al.ll2010l ). Since the tail is rotating, it winds itself into a 
tight spiral, which at larger distances from the pulsar and with limited angular resolution would 
appear as disk-like emission concentrated toward the orbital plane. Moreover, it seems reasonable 
to assume that, in this disk, the pulsar wind becomes well mixed with the dynamically dominant 
companion wind and that both are advected away at the O-star wind terminal speed, v^, at 
dista nces r 3> a (where a is the orbital separation). Adopting Voo = 2400kms _1 ( McSwain et al 



20041 ). we can estimate the time it takes for the wind to reach the angular distance a as td yn 
300(a/l')(uoo/2400 km s" 1 ^ 1 



yr. 



The likely emission mechanism in this case would be synchrotron radiation (see below), im- 
plying the cooling time i syn ~ 300(£'/5 keV) _1 / 2 (l?/10 /iG)~ 3 / 2 yrs. Requiring t sjn > t^ yn places 
an upper limit on the average magnetic field in the extended emission region, 
B < 10(E/5 keV)- 1 / 3 (u oo /2400 km s~ 1 ) 2 / 3 (a/l') -2 / 3 fiG, where E is the observed energy of syn- 
chrotron photons emitted by relativistic electron with Lorentz factor 

7~3x 10 8 (£/5 keV) 1 / 2 ( J B/10/iG)- 1 / 2 . Since the upper limit on B is not much higher than the 
value of the magnetic field in the ISM, it is likely that the actual magnetic field is on the order 
of 10 fiG, and hence the electron Lorentz factor should be of the same order as those in ordinary 
PWNe around isolated puls ars. On the ot her hand, the inferred 7 is close to the maximum possible 
Lorentz factor estimated by iDubusI (120061 ) by requiring the electron gyro-radius to be less than the 
size of the acceleration r egion, a ssume d to be on the order of the standoff distance for the intra- 
binary shock. Moreover, IDubusI (|2006l ) concluded that electrons will not attain 7 > 10 8 because of 
the strong synchrotron radiative losses in the innermost regions. Additional acceleration outside 
the intra-binary shock (e.g., via magnetic reconnection in the spiral-like tail) may provide a way 
to circumvent this limitation. Alternatively, 7 could be somewhat lower if the magnetic field is 
higher than estimated above, which is possible if the bulk flow speed if somewhat higher than the 
terminal O-star wind speed. Indeed, up to a certain distance, the fast pulsar wind could provide 
some re-acceleration to the O-star wind until both components become well mixed. However, the 
detailed modeling of this complex interaction is beyond the scope of this paper. The observed 
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softening of the extended emission with the distance from the point source could be attributed to 
synchrotron cooling which can be substantial at a > 1' if B ~ 10 or somewhat larger. 

Inverse Compton (IC) and adiabatic expansion cooling compete with synchrotron cooling 
and can even dominate, depending on the distance from the particle acceleration region (the 
intra-binary shock). Although at the distances of interest here, the O-star radiation energy 
densit}{§, e ra d = 3.8 x 10~ 10 (a/l') -2 erg cm -3 , would exceed the magnetic field energy density, 
ee = 4.0 x 10~ 12 (i?/10 /^G) 2 erg cm -3 , the IC cooling of 7 > 6 x 10 4 electrons occurs in the 
Klein-Nishina regime and is hence inefficient in comparison with the competing processes (see 
Blumenthal fc Gouldlll97ol l. Indeed, the relevant IC cooling time, nc ~ 20(£/5 keV) 1 / 2 (S/10 /xG)~ 1 / 2 (a/l') 2 Myr, 
on the O-star light will be much larger than the corresponding synchrotron cooling time of electrons 
emitting E = 5 keV photons in B = 10 /j,G. On the other hand, the same electrons will also cool via 
IC upscattering of CMB photons, and this process will dominate over the IC scattering on O-star 
light at the distances of interest. IC scattering off the CMB proceeds in Thompson regime with 
characteristic cooling time ttc,cmb — 7(E/5 keV) -1//2 (-B/10/xG) 1//2 kyr, which is still larger than 
the corresponding r syn (see above). Finally, there is also cooling due to adiabatic expansion, which, 
in 2D geome try, has a ch aracteristic cooling timescale comparable to dynamic timescale (see, e.g. 
Eqn. (28) in bubusl \w0^ ) . 



Zdziarski et al.l (|2010l ) suggested that the X-ray emission from the compact (unresolved) PWN 
in LS I +61 303 could be produced via IC scattering of stellar photons by the pulsar wind electrons 
with 7 ~ 20. The synchrotron emission by the same electrons is expected to be in the radio, while 
the GeV emission could be interpreted as IC up-scattering of stellar light by more energetic photons 
with 7 up to 2 x 10 4 . If the electron SED in LS 5039 extends down to such low 7, these processes may 
contribute significantly at sufficiently small distances from the binary. As the low-energy electrons 
should remain unaffected either by IC or synchrotron cooling on any plausible dynamical timescales, 
one would expect to see radio synchrotron emission on much larger angular scales than it has been 
reported. Deep radio observations, sensitive to extended structures on arcsecond-arcminute scales, 
can provide a useful diagnostics in this case. 



From modeling of absorption and occultation, ISzostek fc Dubud (|201ll ) estimated the plausible 
range of E = (0.1-38) x 10 36 ergs -1 for the alleged pulsar in LS5039. The observed unabsorbed 
luminositj0 Lo.5-8keV = 1 x lO 32 ^ 2 . 5 ergs -1 implies the radiative efficiency range 3 x 10 -6 -10~ 3 , 
which well matches the range f or resolved PWNe around subsonically-moving iso lated Vela-like 



pulsars (lO-MO" 3 ; Fi gure 7 in iKargaltsev et al 



2007 



Kargaltsev fc Pavlovll200cft . If LS 5039 is 



indeed associated with SNR G16.8— 1.1, it has a k inematic age of 40-150 k yr, i.e., a factor of a few 
older than the Vela-like pulsars listed in Table 2 of lKargaltsev et al.1 ([2007). This would favor lower 
values of efficiency from the range given above, and, correspondingly, lower O-star mass loss rates 



following ISzostek fc Dubusl ([201 il l, we assume R+ = 9.3i? and T* = 39, 000 K for the O-star. 

7 Here we take into account only the inner extended emission, because the outer region may include a significant 
dust scattering contribution. 
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from the range discussed in ISzostek Dubusl (|201ll ) 



If the extended emission is due to the pulsar wind, there can be non-variable TeV emission 
associated with the extended wind-filled regio n, even larger than that res olved in X -rays, resembling 



relic P WNe around isolated pulsars (see e.g.. iKargaltsev Pavlovll2010l ). Indeed. I Aharonian et al 



(|2006l ) reported orbital modulation of the LS 5039 TeV flux peaking near inferior conjunction phase 
but also reported a second component with a different spectrum seen near the light curve minimum. 
While the harder emission at the peak of the light curve is likely to come from the innermost region 
of the binary, the soft sp ectrum of the second component is akin to those of TeV PWNe (see 
Kargaltsev fc Pavlovll2010l and references therein) and could be attributed to the extended PWN. 
Therefore, further high-sensitivity, high-resolution observations in TeV are certainly warranted. 



4.3. Transient sources in the field 



We have found several (at least two, see §3.3) transient X-ray sources, all of which are dis- 
tributed within a very narrow strip south of LS 5039. Taken together, the positions of these sources 
at different epochs are inconsistent with a constant velocity or uniformly decelerated motion of a 
single source (see Figures 1 and 7). Therefore, they cannot represent the propagation of a single 
clump of plasma ejected by LS5039. The sources could be unrelated flaring objects (e.g., heav- 
ily absorbed background CVs in the central region of the Galaxy) or, alternatively, all could be 
flare-ups in an outflow largely invisible in X-rays, in which case they don't have to be causally 
connected. Although relativistically outflowing plasma clumps have been seen before in fiQSOs, 
their variability at such large distances from the central source has never been s een as dramatic as 
would appear to be the case here (e.g., iBlundell et al.1 120071 : ICorbel et al.ll2002l ). The direction of 
this outflow would be different from that of the two radio jets seen on milli-arcsecond scales, but 
the latter were seen to change orientation by 12°±3° (jRibo et al.ll2008l ). 



We have searched for other variable sources within the XMM-Newton EPIC-PN FOVs (the 
smallest FOVs among the observations; see Figure 1) and did not find any other variable sources. 
Further monitoring of LS 5039 with Chandra or XMM-Newton is needed to draw any definitive 
conclusions on the nature of the transient sources we discovered. Confirming ejecta clumps would 
favor accretion/ejection scenario over the pulsar wind one. 
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4.4. Comparison with LS I +61 303 and B1259-63 



Among known HMXBs, the one that resembles LS5039 the mosjfl is LSI +61 303. In both 
systems, the donor stars are massive early type stars (O and Be stars with masses 23M Q and 12M Q , 
in LS5039 and LSI +61 303, res pectively) . Both exhibit very similar X-ray properties (including 
fluxes, spectral slopes, Njj; e.g., iRea et al.ll2010l ) and are located at similar distances. Recently, 
both systems were detected in GeV and TeV, where they als o share many common properties 
such as similar VHE spectra and light curves ([Hill et al.l [201CJ, and references therein) . Among 
the differences are a factor of 7 larger orbital period (P rb = 26.5 days) and a factor of 2 larger 
eccentricity (e = 0.72) of the LSI +61 30 3 orbit. Also, LSI +61 303 appears to be more variable in 
X-rays, sometimes showing short flares ([Torres et al.ll2010l ). This could be due to non- uniformity 
of the companion wind which is believed to be concentrated in a rather thin disk in the equatorial 
plane (jSierpowska-Bartosik k, Torres! l2009l ) . The compact object i nteracts with this dense wind 
only during a relatively short phase interval while it crosses the disk (jSierpowska-Bartosik &: Torres 
2009). The wind of the O-star in LS 5039 is b elieved to be much more isotropic. Analyzing a 
49 ks Chandra ACIS observation of LS I +61 303, iParedes et al.1 (120071 ) found faint (58 ± 18 counts) 
asymmetric excess emission extending ~ 12" from the point source. This emission was not confirmed 
in the subsequent longer (96 ks) Chandra ACIS observation, but that observation used Continuous 
Clocking mo de with much hi gher background and only one spatial dimension, greatly complicating 
the analysis (IRea et al.ll2010l ). Finding a much fainter extended emission around LS I +61 303 while 
having the point source X-ray properties very similar to those of the LS 5039 point source, supports 
the conclusion that most of the extended emission in LS 5039 is not a dust scattering halo, if the gas- 
to-dust ratios and dust grain properties are similar for these two systems. If the extended emission 
around the two sources were both dust halos, then, with a similar dust grain size distribution, 
the dust-to-gas ratio would need to be a factor of 10 larger towards LS 5039 in order to explain 
the brighter halo. The dust-to-gas ratio is known to have large scatter along different sight-lines, 
however (see Appendix A). 

If the extended emission around LS 5039 is due to particles supplied by the pulsar wind, one 
could speculate that this emission should be fainter in LS I +61 303 because there the pulsar spends 
most of the time outside the dense companion wind, and hence the pulsar wind is not as efficiently 
confined and decelerated as in LS5039. This means that in LSI +61 303 most of the pulsar wind 
is leaving the binary with large bulk velocity, which should result in much larger but fainter X-ray 
nebula, undetectable in existing data. 

Another object of possibly similar nature is B1259— 63, a HMXB with a much wider orbit 
(-Porb = 3.4 yr), high eccentricity (e = 0.87) and a known pulsar (P = 48 ms, E = 8.3 x 10 35 erg 
s _1 ), where iPavlov et al.l ( 2011al ) found evidence of extended emission. In that case, the emission 



A new HMXB, 1FGL J1018.6— 5 856, has been discov ered very recently (jCorbet et 1201 ll ) whose properties 



are even closer to those of LS 5039 ([Pavlov et al 
LSI +61 303. 



201 lb! ) . but its has not yet been investigated in such detail as 
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is confined to r < 15". This shows that a pulsar system is indeed capable of creating extra-binary 
extended X-ray emission, and therefore that such a scenario may be valid for the case of LS 5039, 
although in this case on a much larger spatial scale. 

5. Summary 

We have discovered asymmetric emission around LS 5039, extending up to 2' from the point 
source. Although it is possible that some of this emission is due to a dust scattering halo, most of 
it is likely produced by energetic particles emanating from LS 5039. Although large-scale energetic 
jets have been observed emanating from ^QSOs, there is a lack of an obvious mechanism to fill a 
large extended volume with X-ray-emitting particles; the pulsar scenario for the compact object 
appears to be preferable. If the wind nebula interpretation of the extended emission is true, one 
would expect to find an even larger radio and/or GeV/TeV nebula, in sensitive enough observations. 
We find transient sources in the field in a narrow stripe south of LS 5039 whose nature is unclear. 

This work is based on data acquired with the Chandra and XMM-Newton X-ray Observatories. 
It was supported by the ACIS Instrument Team contract SV4-74018 (PI G. Garmire) and partly 
supported by archival Chandra grant AR8-9009X, NASA grants NNX09AC84G and NNX09AC81G 
and the National Science Foundation under grants No. 0908733 and 0908611. The work of GGP was 
also partly supported by the Ministry of Education and Science of Russian Federation (Contract 
No. 11.G34.31.0001). 
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Dust halos, often seen around bright point-like X-ray objects, are formed by scattering of 
source X-ray photons on dust grains. Here we will only discuss the case of dust optically thin with 
respect to the photon scattering, T scat < 1, and consider only azimuthally symmetric halos (which 
implies that the dust distribution across the line of sight (LOS) is uniform within the interval of 
angles 6 at which we see the halo). In this case the spectral halo intensity (phcm~ 2 s _1 keV -1 
arcmin -2 ) is given by the equation 



A. 



Dust halo model 




(Al) 



This preprint was prepared with the AAS IATJ5X macros v5.2. 
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where F(E) is the point source spectral flux (photons cm 2 s 1 keV 



(D - d)/D is the 



dimensionless distance from the X-ray source to the scatterer (D and d are the distances from the 
observer to the source and the scatterer, respectively), 9 S ~ 9/x (for small angles) is the scattering 
angle, f(x) is the dimensionless dust density distribution along the LOS (J f{x) dx = 1), and 
da s {E,9 s )/dVL s is the differential scattering cross section per one hydrogen atom, averaged ove r 



Mathis Lee 



1991 



the dust grain distribution over sizes and other grain properties (see, e.g. 
Here we assume that the source spectrum F(E) does not vary: strong variability in the source can 
change the appearance of the halo in a complicated way, as time-lag depends upon radial distance 
from the source. 

To understand the halo properties from simple analytical expressions, we will use the Rayleigh- 
Gans (RG) approximation, in which the total scattering cross section cx E~ 2 ; this approximation 
works better for higher energies, E > 0.5-2 keV, depending on the dust model. For some dust 
models, the averaged differential cross section in the RG can be approximated as (Draine 2003) 



da s {E,9 s ) a s (E) 



<ifL 



< 50 (i + eyei^f 



(A2) 



where 



7 s,50 



S 



10" 22 cm 2 



(A3) 



- and a s (E) « ^ 

are the median scattering angle and the total cross section, respect i vely; E is the energy in keV. 
The constant in first eq. (A3) depends o n the dust m odel; [D raine (200 31) de rived G = 360" from 
the dust mo del of IWeingartner fe Draind (j200ll ). while iBocchino et al.l (120101 ) found = 7.4' for 
the model of ISmith fc Dwekl (|l998l ~ 



It follows from the second eq. (A3) that the scattering optical depth is 

r scat (£) « SN Ht22 E- 2 . 



(A4) 



T he facto r S in the sec ond eq. (A3) is a c o nstan t of the order of 1; e.g., S ~ 1.3 from Figure 6 
of iDraind (|2003l ) , while iPredehl &; Schmittl (|1995l ) found a mean val ue S ~ 0.49 for a number of 



halos observed with ROSAT (but the scatter was very large) , while iMathis &; Led (|199ll ) discuss 
models with S = 0.903, 1.09, and 0.47 (see their Table 1). Costantini (2004; PhD thesis) estimated 
T sca (lkeV) for a number of halo sources observed with Chandra; the values of S derived from her 
results show a very strong scatter, S from 0.018 to 2.26. The scatter itself may be natural, as the 
dust properties may be different for different sources. 

It should be noted that the correlation of r sca (lkeV) with visual extinction Ay. 



r scat (l keV) = (0.056 ± 0.01) A 



(A5) 



(jPredehl h Schmitt 
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1995) is better than that with Njj, but Ay is rarely known for the objects of 
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Substituting (A2) and (A3) in (Al), we obtain the spectral intensity profile 



ar 



. 9E\ 

1 + 1 *e) 



dx . 



(A6) 



For comparison with the point source + halo profile observed in the energy range E\ < E < E 2 , 
the sum of the spectral intensities should be convolved with the detector response, with allowance 
for the image spread caused by the telescope and the detector. We have checked that the energy 
redistribution in the detector only slightly affects the broadband radial profile for a smooth incident 
spectrum. Therefore, assuming the observable halo size to be much larger than the PSF width, we 
obtain 



W0)= dEA cG (E)F(E) 

J E x 



ij(6,E)+N H . 



S f 1 



dx 



1+1 ie) 



(A7) 



where A e g(E) is the detector's effective area, and ip{0,E) is the normalized PSF, which can be 
taken from a simulation (e.g., with MARX for Chandra images). The first term in Equation (A7) 
corresponds to the point source, while the second term describes the halo. 

The above equation can be integrated for a given set of halo parameters, and compared directly 
with the data. In particular, for the dust halo model shown in Figure 4 we picked = 360", S = 1, 
and the dust distribution function 



/(*) 



x < xi, 

(x 2 - Xl)^ 1 Xl < X < X 2 , 
X > X2, 



(A8) 



for xi = 0.1, x 2 = 0.6. According to (A6), this distribution corresponds to the following profile of 
spectral intensity 



'halo 



(e,E) 



F(E)N H>22 S 
2^E\y 2 -y x ) 



arctan 



2/2 -yi (s/2 -yi)(i -ym) 



1 + 2/12/2 (1 + yj)(l + y 2 2 ) \ ' 



(A9) 



where = Xi@/(9E). For small, intermediate, and large 9, Equation (A9) turns into 
F(E)N Hy22 S 



'halo 



(9,E) 



vrB 2 



{xix 2 ) 

(tt/4x 2 )(@/9E) 



9 < Xl e/E, 
Xi q/e < 9 < x 2 q/e, 



(A10) 



(l/3)(e/9E) 4 (xj + Xl x 2 + x 2 2 ) 9 » x 2 e/E, 



where the approximation for the intermediate 9 implies x 2 S> x\. It follows from Equation (A10) 
that, for the model (A8) (which may correspond to the case when the scattering occurs mainly in a 
Galactic arm), the halo profile consists of three parts: a flat top (whose size is proportional to xi), 
a slowly decreasing part (oc 9~ r ), and a steeply decreasing part (oc 9~ A ), so that the characteristic 
size of a halo is about x 2 Q/E. 



